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The analysis of a physical problem simplifies considerably when one uses a suitable coordinate 
system. We apply this approach to the discrete-time quantum walks with coins given by 2j -|- 1- 
dimensional Wigner rotation matrices (Wigner walks), a model which was introduced in T. Miyazaki 
et al., Phys. Rev. A 76, 012332 (2007). First, we show that from the three parameters of the coin 
operator only one is physically relevant for the limit density of the Wigner walk. Next, we construct 
a suitable basis of the coin space in which the limit density of the Wigner walk acquires a much 
simpler form. This allows us to identify various dynamical regimes which are otherwise hidden in 
the standard basis description. As an example, we show that it is possible to Hnd an initial state 
which reduces the number of peaks in the probability distribution from generic 2j -|-1 to a single one. 
Moreover, the models with integer j lead to the trapping effect. The derived formula for the trapping 
probability reveals that it can be highly asymmetric and it deviates from purely exponential decay. 
Explicit results are given up to the dimension hve. 


I. INTRODUCTION 

Quantum walk m is an extension of the concept of a 
random walk ^ to unitary evolution of a quantum par¬ 
ticle on a graph. Both continuous-time and discrete-time 
(or coined) quantum walks have found promising appli¬ 
cations in quantum information processing, e.g. in quan¬ 
tum search algorithms EHH], graph isomorphism test¬ 
ing [MTT] , finding structural anomalies in graphs p'2HT4] 
or as universal tools for quantum computation [13 US]. 
Quantum walks were realized in a number of experiments 
using optically trapped atoms m, cold ions (TS] |T9| and 
photons [2nH26| . 

The dynamics of the discrete-time quantum walks cru¬ 
cially depends on the choice of the coin operator and 
the initial state which has essential implications for 
the performance of the search algorithm |28|. The be¬ 
haviour of two-state quantum walks with homogeneous 
coin is fully understood [5S| . This is not true for quantum 
walks with higher-dimensional coins, since the complex¬ 
ity of the unitary group grows rapidly with its dimen¬ 
sion. Most of the known results were obtained for quan¬ 
tum walks with Grover coin which leads to the so-called 
trapping (or localization) effect [33 El] ■ In contrast to 
the two-state walk, which show purely ballistic spread¬ 
ing, the particle performing the Grover walk has a non¬ 
vanishing probability to remain trapped in the vicinity 
of the origin. This is represented by a stationary peak in 
the probability distribution located at the origin, which 
does not vanish with the increasing number of steps but 
decays exponentially with the distance from the starting 
point. The Grover walk and its extensions were inten¬ 
sively studied, either for line [SOUHH] - plane [37U89] or 
higher dimensional lattices [13 SI]- 

Among the few models of quantum walks on a line 
which were solved analytically for arbitrary dimension 
of the coin operator, i.e. arbitrary number of displace¬ 
ments particle can make in a single time-step, are the 
Wigner walks introduced in [12]. This model utilizes 


2j + 1-dimensional Wigner rotation matrices as coin op¬ 
erators and is closely related to the tensor product model 
of quantum walks |431146j . Explicitly, using the Fourier 
analysis m and the weak-limit theorem [48j . the limit 
density of the Wigner walks was derived in [12] for arbi¬ 
trary j. However, as the authors point out in [42] . the 
dependence of the limit density on the initial coin state 
is rather involved. The reason for this inconvenience is 
that the initial coin state is decomposed into the stan¬ 
dard basis of the coin space, i.e. the one in which the step 
operator is defined. The aim of our paper is to present 
an alternative approach to |42| . Namely, by using a more 
suitable basis of the coin space we put the results of [12] 
into much more convenient form. This allows us to iden¬ 
tify previously unknown features which are hidden in the 
standard basis description. As an example, we show that 
while the generic probability distribution of the 2j + 1- 
dimensional Wigner walk has 2j + 1 peaks, there exist 
initial conditions for which the number of peaks reduces 
to one. Moreover, the models with integer j exhibit the 
trapping effect similar to the Grover walk. This feature 
was not analyzed in detail in [12]. We derive the ex¬ 
plicit formula for the trapping probability for j = 1 and 
j = 2. Our results show that the trapping probability 
can be highly asymmetric. Moreover, the decay of the 
trapping probability with the distance from the origin is 
not exactly exponential for j = 2. 

The paper is organized as follows. In Section [IT] we 
review the results on quantum walks with Wigner coin 
obtained in [42] . We find that when we focus only on the 
limit density and its moments the number of physically 
relevant coin parameters can be reduced from three to 
one. Next, we outline the procedure to determine the 
suitable basis of the coin space in which the limit den¬ 
sity simplifies considerably. Namely, we choose the basis 
states among those which give rise to non-generic distri¬ 
butions. In Section jin] the suitable basis is determined 
for the simplest case of the two-state model. Moreover, 
we find a recipe for the construction of suitable bases 
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for higher-dimensional models. We apply this recipe in 
Section ITVI to the three-state model and we find that it 
is equivalent to the three-state Grover walk [301131] and 
its one-parameter generalization |53H55] . In Sections [V| 
and |VI| the four-state and five-state models are studied, 
respectively. Expressing the limit density in the suitable 
basis allows us to determine additional non-generic situ¬ 
ations. In particular, we find initial conditions such that 
the number of peaks in the probability distribution is 
reduced to one. Finally, we summarize our results and 
present an outlook in Section [Vll[ More technical details 
of the analysis are left for the Appendices. 


II. QUANTUM WALKS WITH WIGNER COINS 

In this Section we review the quantum walks on a line 
with Wigner coins (Wigner walks) following their intro¬ 
duction and analysis in [42]. The Hilbert space of the 
Wigner walk is given by a tensor product 

n = np^nc, 


and the coin space He- The dimension of the coin space 
is 2j + 1, where j is a (half-)integer. The standard basis 
of the coin space is formed by vectors \m) corresponding 
to the jumps of length 2m where m = — j, — j + 1 ..., j. 
Single step of the time-evolution is given by a unitary 
operator 

17 = 5-(JGi?«(a,/3,7)), (1) 

where S is the step operator 
-1-00 j 

E E |a;-f 2 m)(a;| G |TO)(m|, ( 2 ) 

X——00 m——j 

and (a, (3, 7 ) denotes the coin operator which is given 
by the Wigner rotation matrix [49] [50], i.e. the 2j + 
1 -dimensional irreducible representation of the rotation 
group S0(3). The matrix elements of the coin in the 
standard basis of the coin space 

■Rmn(a>/3>7) = (w|.R('^^(a,/3,7)|n), 

are given by 

rW(/3)e-*^”, 


of the position space 

Tip = Span{|a;)|a; S Z} . 


= e“ 


where 


r-0) (^) = ^r(j,TO,n,0 ("cos I") 


2j-\-m—n—2l ^ 2/ —m+n 


V 2 / 


sm 


The factor T(j, m, n, 1) reads 


T{j,m,n,l) = (-1)' 


V{j + my.jj - m)\{j + n)\{j - ri)! 

{j — n — l)\{j + m — l)l{l — m + n)\l\' 


(3) 


The summation index I runs over all integers such that 
all factorials in ^ are well defined. 

As the initial state of the Wigner walk we consider 

|iA(0)) = |0) G IV’c), 

i.e. the particle starts from the origin with the initial 
coin state 

j 

\^c) = E 9 m\m). (4) 

The amplitudes qm fulfill the normalization condition 

E = 1- 


The state of the particle after t steps of the walk reads 

which can be decomposed into the standard basis as 

IV’W)=E E (^\m). 

X m——j 

The probability to find the particle at position x after t 
steps of the quantum walk is then given by 

m— — 

We note that all Wigner walks are bipartite, i.e. half of 
the lattice points are empty at any time step. For integer 
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j the odd sites are never occupied, while for half-integer 
j the walk oscillates between even and odd sites. 

Since the Wigner walks are translationally invariant, 
their analysis is greatly simplified in the Fourier picture 
m- Moreover, in the asymptotic limit the moments of 
the particle’s re-scaled position (or pseudo-velocity) can 
be expressed in the form [48] 

where {v) is the limit density. Its explicit form for 
Wigner walks was derived in |42j where it was shown 
that it is given by a sum 

u^i\v)= ( 6 ) 

0<m<j 

The summation index m runs over (half-)integers, de¬ 
pending on whether j is half-integer or integer, in unit 
steps. The individual densities have the form 



where n{v,a) is the Konno’s density function [5T1152] 


n -VT-r’ \ ! 

7r(l — v^)y'[a — v)(a + v) 

The symbol Af (f) denotes the weight function which 

is a polynomial of degree 2j in v with coefficients deter¬ 
mined by the initial state and the coin parameters j3 and 
7 . The decomposition W of the limit density 
into [_)] densities (uj indicates that the 2 j-|- 1 -state 
Wigner walk can be considered as superposition of |"j] 
walks which propagates through the lattice with differ¬ 
ent velocity given by 2m cos ^. In addition to the evalu¬ 
ation of moments the limit density §) can be used 
to approximate the shape of the probability distribution 
of the Wigner walk in finite time according to 

(9) 

The factor of 2 accounts for the fact that the Wigner 
walks are bipartite. Note that the Konno’s density func¬ 
tion ^{y]a) diverges for v —?> ±a. These divergencies 
correspond to the 2 j - 1 - 1 peaks ( 2 j if j is an integer) in 
the probability distribution p^^\x,t) which are found at 
X « ± 2 m cos where the range of the index m is the 
same as in the sum Q. 

We note that for integer j the particle is allowed to 
stay at its actual position, while for half-integer j it has 
to leave the previously occupied site. This has a crucial 
impact on the dynamics of the quantum walk, as was 
found in |42|. Indeed, for half-integer j the spectrum of 
the evolution operator 0 is purely continuous, while for 
integer j it has, in addition, an isolated eigenvalue with 


infinitely-many localized eigenstates. The presence of the 
point spectrum results in an additional peak in the center 
of the probability distribution which decays rapidly with 
the distance from the origin. Moreover, the peak does 
not vanish with increasing number of time steps, i.e. for 
integer j the particle has a non-vanishing probability to 
remain at position x in the asymptotic limit. This feature 
was first observed for three-state Grover walk on a line 
[sniEi]- We denote the limiting value of the probability 
to remain at position x as 

= \im p^^\x,t), ( 10 ) 

and call it the trapping probability. Its explicit form was 
not given in [42]. We evaluate the trapping probability 
in Sections jW] and lYg where we treat the Wigner walks 
with j = 1 , resp. j = 2 , following the approach used in 

EMSSj. 

In principle, the formulas 0, 0 and the general ex¬ 
pression for the weight functions derived in 

[35] allow us to calculate the limit density for ar¬ 

bitrary (half-)integer j. However, the actual form of the 
limit density is rather complicated already for small val¬ 
ues of j, as can be observed in [35] where the authors 
provided the explicit results up to j = |. This makes 
further investigation of Wigner walks quite demanding. 
It is the aim of the paper to simplify the limit density as 
much as possible. In order to do so, we first discuss which 
of the coin parameters are physically relevant. Second, 
we simplify the dependence of the weight function on the 
initial coin state by choosing a more suitable basis in the 
coin space. Since we are interested only in the moments 
and the limit density, we consider two Wigner walks with 
coins .R^'^^(q;i,/ 3 i, 7 i) and {a 2 , (32,72) as equivalent if 
for any initial coin state |' 0 Ci) of tho first walk there ex¬ 
ists an initial coin state of the second walk \i(’C 2 ) such 
that the resulting limit densities are the same. We note 
that two equivalent Wigner walks, as we have defined 
them, might exhibit different properties when additional 
aspects of quantum walks, such as topological phases 
l60| . are of interest. 

Let us begin with the coin parameters. Notice that 
Konno’s density function depends only on one of them, 
namely (3. Moreover, as was shown in [31] the weight 
functions Af are determined by the initial coin state 
and the coin parameters (3 and 7 , but they are indepen¬ 
dent of a. Hence, all quantum walks with different a are 
equivalent, as far as the limit density and the moments 
of the particle’s position are concerned. Therefore, we 
consider a = 0 from now on. Moreover, the dependence 
of the weights on 7 is rather simple. As was 

found in |42j the parameter 7 enters only through the 
terms of the form Hence, the Wigner 

walk with coin (0, (3, 0) and the initial coin state 0 
gives the same limit density as the Wigner walk with the 
coin ( 0 ^ ^) and the initial coin state 

\^c) = H 

m^-j 
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Therefore, all models with different 7 are equivalent 
according to our definition and we restrict our fur¬ 
ther analysis to the case 7 = 0 . Hence, we are left 
with only one physically relevant coin parameter /3. 
We note that for two-state walks the equivalence of 
the three-parameter set of quantum walks with coins 
^( 1 / 2 ) (q,^ 7), which in fact covers all translationally 

invariant two-state quantum walks, with the single¬ 
parameter family /3, 0) was already established 

in [6l]. The above discussion shows that within the set 
of Wigner walks the equivalence hold for arbitrary j. 

We now turn to a slightly different parametrization of 
the coin operator /3, 0) which is more suitable in 

the context of Wigner walks. Namely, instead of /3 we 
consider the parameter p given by 

P 

P = cos -, 

which corresponds to the divergencies of the limit densi¬ 
ties Q. It also directly determines the speed of propa¬ 
gation of the wave-packet through the lattice |S2]. Since 
the Euler angle /3 is limited to the interval [0, tt], the new 
parameter p varies from 0 to 1, and the identity 

sin ^ 

holds. Finally, we define the coin operator R^^\p) = 
j5,Q) with the matrix elements in the standard 
basis given by 

R^nip) = {m\U^\p)\n) =^r(j,TO,n,Z) 

I 

/ , _\ 2l — m-\-n 

^^2j+m-n-2l . (11) 

Let us now turn to the dependence of the limit density 
on the initial coin state. As was pointed out in [42], the 
weight functions are rather involved functions of 

the coefficients qi of the initial coin state in the standard 
basis. The aim of this paper is to simplify these expres¬ 
sions by choosing a more suitable basis of the coin space. 
Following the common experience of quantum mechanics 
one would expect that the suitable basis is the one given 
by the eigenvectors of the operator involved. This is in¬ 
deed the case of the three-state Grover walk |35|. How¬ 
ever, for Wigner walks the eigenvectors do not represent 
the best choice. To construct the suitable basis we con¬ 
sider one additional property of the eigenvectors of the 
Grover coin used in |35j . Namely, the basis vectors were 
chosen such that they result in non-generic distributions 
of the Grover walk, i.e. they reduce the number of peaks 
in the probability distribution. We adopt this require¬ 
ment for Wigner walks and select the suitable basis vec¬ 
tors among such states. The conditions for non-generic 
distributions are straightforward in the asymptotic limit, 
since then the peaks correspond to the divergencies of 
the limit density (|^. Therefore, we have to determine 
the states for which some of the divergencies of Q van¬ 
ish. In the following Section we solve these conditions 


and determine the suitable basis states for the two-state 
model. We then rewrite the suitable basis in terms of the 
eigenvectors of the coin operator. This gives us a recipe 
for construction of suitable bases in higher-dimensional 
models. 


III. TWO-STATE MODEL 


We begin with the simplest case of a two-state model 
where J = |. The coin space is two-dimensional with 
the basis states |l/2), | — 1/2), which correspond to the 
jumps of one step to the right and one step to the left. In 
the standard basis the matrix representation of the coin 
operator is given by 

The limiting probability density calculated in [35] reads 

where p{v; p) is the Konno’s density function Q and the 
weight is given by 


A1 (i/ 2’1/2) = l-f Alli/"’i/')u. 


(14) 


(■ 2^ /2 2 /2) 

The linear term A4) in the weight has the form 


^(1/2,1/2) ^ + |(?-l/2r + 


.yiE 

p 


(91/29-1/2 + 91/29-1/2)1(15) 


where g/s represent the coefficients of the initial coin 
state in the standard basis of the coin space 


li ^ c ) — 91/2I1/2) + 9-1/2I —1/2). 


(16) 


We will now determine the suitable basis for the descrip¬ 
tion of the two-state model. Following the discussion in 
the previous section, we will construct the suitable basis 
from the states which lead to non-generic probability dis¬ 
tributions, i.e. the states for which one of the divergen¬ 
cies of the limit density (13) vanishes. These divergencies 
coincide with those of the Konno‘s density function (§, 
which appear for v = ±p. To eliminate them we have 
to find such q±i /2 that the weight function tends to zero 
faster than the denominator of the Konno’s density func¬ 
tion for V = ±p, i.e. the weight has to attain the form 


Al(5’5) = l±-. 

P 
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Hence, the linear term (15) has to be equal to ±i. The 
solutions of these two equations 


1-^ 


1 + - 
P 


<Zl/2 


9-1/2 


91/2 


9-1/2 — 


1+ P 


1 - p 


I- p 


1+ P 


provide the coefficients of the initial states (161 that elim¬ 
inate the divergence of the limit density (13) at i; = ±p. 
We denote these vectors as |x^)- In the standard basis 
they have the form 


lx+) = 


lx ) = 


\ + p 


I-p 


|l/2) - 


|l/2) + 


P\ 


1 + p. 


- 1 / 2 ), 


- 1/2). (17) 


Clearly, these two states form an orthonormal basis of 
the coin space. Moreover, we find that in this basis the 
weight function (14) simplifies considerably. We express 
the initial coin state in terms of the basis {|x^)} accord¬ 
ing to 


IV'c) = h+\x") + h lx ). 


The correspondence between the amplitudes of the initial 
coin state in the suitable basis hi and in the standard 
basis is then given by 


9i/2 


h-l/2 


1 + P 


h+ + 


2 


I- p 


h+ 


2 


Inserting these relations into the formula (15) we find 


that the linear term of the weight function reduces into 

Hence, the limit dens ity for a two-state quantum walk 
with Wigner coin (12) in the suitable basis {|x^)} reads 




^r^(i + (i-2|/i+p)g) 

7r(l - v'^)yj{p - V){p + V) 


(18) 


This result shows another benefit of using the basis 
{|x^)}. Namely, by absorbing part of the p dependence 
into the definition of the basis states ((Tt]), p is effectively 
reduced to a scaling parameter. Indeed, varying the value 
of p while keeping the amplitudes and h~ untouched 
does not affect the shape of the limit density (18). The 


parameter p simply determines how far the density is 


stretched. We point out that the same will apply for 
Wigner walks with higher values of j. 

To illustrate our results we display in FIG. the prob¬ 
ability distribution for the initial coin state |x~) after 
100 steps. The red curve depicts the limit density which 
for |x“) is given by 


V 


( 1 / 2 ) 

lx-> 


(u) 


p^y/p + v 
7rp(l — v'^)y/p — V 


Clearly, the limit density diverges for u —)■ p but tends 
to zero for v —> —p. This corresponds to the presence of 
only one peak in FIG. which we emphasize by using a 
logarithmic scale on the y-axis 0 



X 


FIG. 1. (Color online) Probability distribution after 100 steps 
of the two-state Wigner walk with the coin (p). The ini¬ 

tial state was chosen as |x~). The coin parameter is p = 0.8. 
We note that for different values of p the probability distri¬ 
bution will have the same shape, only its width will change 
accordingly. The same applies to all figures. 


Let us point out that the limit density of any homo¬ 
geneous two-state quantum walk on a line can be trans¬ 
formed into the form (18). Indeed, the two-dimensional 
Wigner coins /3,j) covers all SU(2) matrices 

and the equivalence of a general SU(2) walk to the walk 
with the coin was already pointed out in [61| . 

Hence, for any two-state quantum walk on a line one can 
find two orthogonal coin states giving rise to non-generic 
distributions with only one peak. 

Finally, let us rewrite the suitable basis {|x^)} terms 
of the eigenvectors of the coin operator (12). This decom¬ 
position will be useful in the following Sections where we 
treat higher-dimensional models. The eigenvectors of the 
coin operator (12) are given by 

|V>±) = ^(|1/2)t *1-1/2)). (19) 


^ We use log-scale in all figures to emphasize the peaks and dips. 
The black points are given by the numerical simulation of the 
quantum walk. The full red curves always correspond to the ap¬ 
proximation of the probability distribution with the limit density 
given by 
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They satisfy the eigenvalue equations 




where the phase of the eigenvalues reads 


eigenvectors (191 is easily found to be 


lx+) = ^(e-*"|V'+)+e*"|V'-)), 
lx”) = ^ (e”*" 1 ^/’+) - e* 2 1 ^-)^ . 


( 20 ) 


These relations provide us with a recipe for construction 
of suitable bases for j > 


ip = arccos p. 


IV. THREE-STATE MODEL 


The decomposition of the suitable basis {|x^)} into the 


In this Section we analyze the three-state Wigner walk 
which is driven by the coin operator 




_ -V^pyji - p^ 

\/2pa/ 1 - _i + 2p^ 

1 - yj2p^jl - 


l-p2 \ 

/ 


The Wigner matrix (p) is reminiscent of the modified 
Grover coin which was introduced in |32j and recently 
analyzed in more detail in [34] [35] . Therefore, it is not 
surprising that the two models yield the same results, 
as we identify in the following subsections. In fact, the 
affinity of the three-state Wigner walk and the three-state 
Grover walk was already discussed in [42] . 


In the three-state Wigner walk the particle is allowed 
to remain at its present position. As we have discussed 
before, this implies that the probability distribution of 
the three-state Wigner walk is not described solely by the 
limit density ([^. There is an additional non-vanishing 
and stationary peak at the origin due to the trapping ef¬ 


fect (10). We analyze the limit density in subsection IV A 
and the trapping probability in subsection jlVBj Before 
we give their explicit forms, we first construct the suitable 
basis of the coin space. We adopt the recipe provided by 
Eq. (20), which gives the decomposition of the suitable 
basis states for the two-state model into the eigenvectors 
of the coin operator. For the three-state model, the coin 
R^^\p) has two eigenvectors satisfying the eigen¬ 

value equations 


i?W(p)|V'±) = e±*^|V'±). 


The explicit form of the eigenvectors and the phase <p 
is left for the Appendix Following the formula (20) 
for the two-state model, we construct from two or¬ 
thonormal vectors |x^) which will serve as part of the 
new basis. Moreover, the coin R^^^p) has an addi¬ 
tional eigenvector |^o) corresponding to the eigenvalue 
one. This vector is indeed orthogonal to both and 
Ix^)- Hence, we consider it as the last vector of the new 


basis which reads 

Ixo) = IV'o), 

lx+) = ^(e-“"|V'+)+e*"lV'-)), 

lx”) = (e"*2 |^+). (21) 

The initial state is now decomposed in the suitable basis 
{|Xo), Ix"^)} according to 

IV’c) =^o|xo) + /i+lx^)+^"lx")- (22) 

The explicit correspondence between the amplitudes of 
the initial state in the suitable basis hi and in the stan¬ 
dard basis qi is given in the Appendix[^ In the following 
we show that the limit density and the trapping probabil¬ 
ity obtain much more convenient forms when expressed 
in the basis {|xo), Ix"^)}- 


A. Limit density 


The limit density for the three-state Wigner walk was 
derived in [35] and reads 

where can be expressed as a polynomial of de¬ 

gree two in V 


The individual terms depend on the coin param¬ 

eter p and the initial coin state. Their explicit form in 
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the standard basis was given in |42] . we present it in the 
Appendix^for comparison. There we also show how the 
terms M.\ ' ^ simplify when we expand the initial state 
into the suitable basis according to (22). Indeed, in the 
suitable basis we obtain 


= \h+\^ + \h-\\ 

+ hoh~), 


For illustration we show in FIGs. |2|3| and the prob¬ 
ability distribution of the three-state Wigner walk with 
the initial coin state given by one of the suitable basis 
vectors (21). In FIG. we find that for |xo) the limit 
density tends to zero at the origin. Indeed, for ho = 1 
the constant term of the limit density vanishes. 

Nevertheless, the probability distribution at the origin 
does not vanish due to the trapping effed]^ For |x+) both 
peaks at —)■ ±2p vanishes, as we illustrate in FIG. [3 Fi¬ 
nally, FIG. [^indicates that for the last basis vector x~) 
the trapping at the origin disappears, which we identify 
analytically in the Appendix [A| 



X 


FIG. 2. Probability distribution after 100 steps of the three- 
state Wigner walk. The initial coin state was chosen as |xo) 
and the coin parameter is p = 0.5 


We point out that the limit density of the three-state 
Grover walk and its one-parameter extension exhibit all 
of these features [5S]. Moreover, in the following subsec¬ 
tion we show that also the trapping effect of the three- 
state Wigner walk is the same as for the three-state 
Grover walk. 


^ We calculate the explicit form of the trapping probability in the 
following subsection. The trapping probability is in this and all 
relevant figures denoted by the dashed blue curve. 
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FIG. 3. Probability distribution after 100 steps of the three- 
state Wigner walk. The initial coin state was chosen as |x^) 
and the coin parameter is p = 0.8 
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FIG. 4. Probability distribution after 100 steps of the three- 
state Wigner walk. The initial coin state was chosen as |x~) 
and the coin parameter is p = 0.6 


B. Trapping probability 


Let us now turn to the trapping probability. The de¬ 
tails of the calculations are left for the Appendix Fol¬ 
lowing the same approach that was used for the three- 
state Grover walk in [3SH3S] we find that in the suitable 


basis (21) the trapping probability reads 

Q 2 |d?(I^|/jo_ft+| 2 ^ ir<0, 

^((l-p2)|/io|2 +|/j+|2)^ 3-^0, (23) 


Here we have denoted 




2 - p2 - 2^1 - p2 


(24) 


We note that the trapping probability for the one- 
parameter extension of the three-state Grover walk has 


























exactly the same form m- 

In the suitable basis (21), the trapping probability (23) 
is independent of the amplitude h~ of the initial coin 
state (22). This fact reduces the dependence of the trap¬ 
ping probability to just two amplitudes Hq and . How¬ 
ever, the expression (23) can be simplified further by an 
additional change of basis. Notice that the dependence of 
the trapping probability is different for positive and neg¬ 
ative positions. In fact, one can choose such an initial 
state that the trapping effect appears only on positive 
or negative semi-axis. We note that the same feature 
was identified for the three-state Grover walk [33] and its 
one-parameter extension |3S] . For the three-state Wigner 
walk with the initial coin state 
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we find that the trapping probability turns into 


q 2N2(1^|^_|2^ cr< 0 , 

Poo( 2 a;)=< - 2^1^+, x = 0 , 


x>0. 


The advantage of the basis {|x~)j 1^^)} is that the trap¬ 
ping probability outside the origin has a simpler form. 
Namely, it depends only on one amplitude (resp. l~) 
for positive x (resp. negative x). This additional change 
of basis becomes crucial in Section lvTl where we treat the 
five-state Wigner walk. 


|A+) = ^(lxo) + lx+)), 


V. FOUR-STATE MODEL 


the trapping effect appears only for non-negative posi¬ 
tions. Similarly, the state 

|A-) = ^(|Xo)-|x+)) 

shows trapping only for non-positive x. We illustrate this 
feature in FIG. El 


10- 




a, 


10-^ 




\i>c) = lA > 
p = 0.5 


-10 -5 0 5 10 


X 


FIG. 5. Probability distribution after 10000 steps of the three- 
state Wigner walk. The initial coin state was chosen as |A“) 
and the coin parameter is p = 0.5. Only a small vicinity of 
the origin is displayed to emphasize that for | A“) the trapping 
effect appears only for x < 0. Higher number of steps in 
comparison with other figures was chosen so that the trapping 
effect is sufficiently pronounced further from the origin. 


The vectors |A+) and |A“) are mutually orthogonal. 
Moreover, they are both orthogonal to \x~) and the 
triplet {|x~)!|A^)} forms an orthonormal basis of the 
coin space. When we decompose the initial coin state as 

\i’c) = h lx ) -I- Z''"|A'*') -I- Z |A ), 


For J = I we obtain a four state quantum walk model 
with the coin operator determined by 0- The coin 
operator has two pairs of eigenvectors ) and 

with conjugated pairs of eigenvalues 

The explicit form of the eigenvectors and the phases (£]9 
is given in the Appendix]^ Following the formula (20) 
we construct the suitable basis by combining only the 
eigenvectors corresponding to the same phase factor (pi 
or Lp 2 - Therefore, we consider the suitable basis in the 
form 

1 

3' 






IxD = 
Id) - 
I.P = 




'IV’i ) - e 
■|V'2+) + 
\^t)- 


IV'l )), 

e*'^IV'2-))- 


(25) 


The initial coin state is decomposed in the suitable basis 
according to 

He) = ht\xt) + K\xi) + ht\xt) + Ix^)- 


We present the explicit relation between the amplitudes 
in the suitable basis hi and the standard basis Qi in the 
Appendix [^ 

Let us now turn to the limit density which for the four- 
state model is given by a sum of two densities [42] 

l/(3/2)(x) = zz(3/2.1/2)(^) +^(3/2,3/2)(^)^ 


The individual densities corresponding to a slower walk 
(to = 1) and a faster walk (to = |) read 

,.<3/2.3/2)(„) ^ Ip ^(3/2,3/2) Q) _ 
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The weight functions (f) are given by cubic 

polynomials in v 


7W(3/2,m)(^) ^ (27) 

fe =0 


with coefficients determined by the initial coin 

state and the coin parameter p. Their explicit forms in 
the standard basis of the coin space are given in [12] . We 
express them in terms of the suitable basi s {|x f), 1x^1} 
in the Appendix]^ Using the expressions ( |B3| ) and (B4) 
one can show by direct computation that each of the 
basis states eliminates two peaks of the probability 
distribution - one in each of the individual limit densities 
j^(3/2,i/2)(^^ and ^(3/2,3/2)(^^^ rpp illustrate our results 
we present in FIGs. [^andj^the probability distributions 
for the initial states lx)*") and \xt)- 




FIG. 7. (Color online) Probability distribution after 100 steps 
of the four-state Wigner walk. The initial coin state was cho¬ 
sen as \xt)- The coin parameter is p = 0.5. For the initial 
state 1 x 2 ^) thfi probability distribution is a mirror image of 
the present one. 


which lead us to the relations 



-I- 

1 I 



(29) 


Hence, we find a two-dimensional subset of states 


= ht\xt)+K IXi ) 



73 


1x1), (30) 


X 


FIG. 6. (Color online) Probability distribution after 100 steps 
of the four-state Wigner walk. The initial coin state was cho¬ 
sen as IXi")- The coin parameter is p = 0.8. The initial state 
Ixr) results in a distribution which is a mirror image of the 
present one. 


We point out that it is possible to eliminate both peaks 
for each limit density i/G/z.m) (^^ order to 

do so we have to find an initial coin state such that both 
divergences provided by the Konno’s density p{y\p) in 
(26) vanishes. Hence, the weight function have to be of 
the form 


for some arbitrary constants a and b. This is satisfied 
provided that 

_^(3/2.m) ^ p2^{3/2,m) ^ 

_^(3/2.m) p2_y^G/2.m) ^ ^ 28 ) 

Adding and subtracting these two equations for the 
weight function Ad^3/2.i/2)(.^^ obtain the following 

\h++ V^h^\‘^ =0, |/iC= 0, 


for which the peaks of the slower walk described by 
t^(3/2,i/2)(t, ) va nishes. We note that for the vectors of 
the family (30) with either = 0 oi = 0 

in addition one of the peaks of the faster walk vanishes, 
i.e. these states result in a probability distribution with 
only one peak. We illustrate this result in FIG.[^ 



X 


FIG. 8. (Color online) Probability distribution after 100 
steps of the four-state Wigner walk. The initial co in s tate 
-^|xi~) ~ llx^) "'S'S chosen as a part of the set (301 for 
which the inner peaks in the probability distribution vanishes. 
Moreover, since hf = = 0, only one of the outer peaks is 

present. The coin parameter is p = 0.6 
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Similarly, we can cancel both divergencies for the sec¬ 
ond limit density corresponding to the faster 

pea ks in the probability distribution. From the equations 
(281 for m = I we find that both faster peaks vanish pro¬ 
vided that 


h+ = VSh^, h-=V3h^. (31) 

Hence, we obtain a two-dimensional set of initial states 

+ ^rixr) + + v ^ k \ x 2 ), 

(32) 

for which the outer peaks in the probability distribution 
vanishes. Moreover, the vectors from the family ( |3^ sat¬ 
isfying either hf — = 0 or = 0 lead to 

elimination of one additional peak of the slower walk, i.e. 
they result in a probability distribution with only one 
peak. We illustrate this result in FIG.[^ 



X 


FIG. 9. (Color online) Probability distribution after 100 steps 
of the four-state Wigner walk. The initial coin state \ \xt) + 
IxJ) is an element of the set (321 for which the outer peaks 


in the probability distribution vanishes. Since = hj = 0, 
only one of the inner peaks is present. The coin parameter is 
p = 0.8 


Clearly, it is impossible to fulfill both conditions (29) 


and (31) simultaneously. Hence, the probability distri¬ 


bution of the four-state Wigner walk has always at least 
one peak. 


of conjugated eigenvalues 


The explicit form of the eigenvectors and the phases Lpi ^2 
is left for the Appendix From these two pairs of eigen¬ 
vectors we construct four basis vectors |xf) and |xj) 
according to (25). In addition, the coin has an 

eigenvalue 1. We include the corresponding eigenvector 
lipo) as the last basis vector. Hence, we obtain the suit¬ 
able basis in the form 


Ixo) = IV'o), 

\xt) = Te^^lV^r)) > 

Ixr) = ^ , 

Ix^) = ^ , 

1X2-) = ^ (e-*^|^2+) • (33) 

The initial state is decomposed into the suitable basis 
according to 

IV'c) = ^olxo) + h+\xt) + K\xl) + htlxt) + Ix^)- 


We present the correspondence between the amplitudes 
in the suitable basis hi and in the standard basis qi in 
the Appendix In the following subsections we illus¬ 
trate that the change of the basis allows us to identify 
various interesting regimes of dynamics which are other¬ 
wise hidden in the standard basis description. 


A. Limit density 

Let us first discuss the limit density. As for the four- 
state walk of Section |vj the total limit density ([^ is a 
sum of two densities 

corresponding to a slower walk 


VI. FIVE-STATE MODEL 2 V2 / V2/ 


In this Section we analyze the case j = 2 which leads 
to the five-state Wigner walk. Similarly to the three- 
state model, which we have discussed in Section [TV] the 
probability distribution of the five-state model consists 
of the limit density and the trapping probability. We 
treat them in subsections VIA and |VIB[ but first we 
construct the suitable basis of the coin space following 
the recipe ( [2o| . The coin operator has two pairs 

of eigenvectors and IV'^) corresponding to two pairs 


and a faster walk 

The weight functions AI*'^’’”^(u) are polynomials of de¬ 
gree four in v 
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with coefhcients determined by the initial coin 

state and the coin parameter p. Their explicit form in the 
standard basis of the coin space can be evaluated using 
the procedure of [55]. In the Appendix we express 
them in terms of the coefficients in the suitable basis. 

Let us now illustrate the role of individual vectors of 
the suitable basis {|xo)) Ixf )> Ix^)} oii the dynamics of 
the five-state Wigner walk. Using the explicit form of 
the weight functions (C3) and (C4) one can show that 
each of the basis states eliminates two peaks, either in 
the slower walk described or in the faster walk 

given by In FIG.j^we display the probability 

distribution for the initial state |xo)- We find that both 
peaks of the slower walk vanishes. In addition, both den¬ 
sities and tend to zero at the origin, 

since both weight functions miss the constant 

term In FIG. this corresponds to the signif¬ 

icant dip of the red curve around the origin. However, 
the total probability distribution does not vanish at the 
origin due to the trapping effect, which is illustrated by 
the blue curve in FIG. [101 



X 


FIG. 10. (Color online) Probability distribution after 100 
steps of the five-state Wigner walk. The initial coin state was 
chosen as |xo). The coin parameter is p = 0.4 


In FIG. 
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we have chosen the initial state as |x^)- For 


this particular state both outer peaks vanishes. FIG. 12 


displays the probability distribution for the initial state 
Ixi”)- In. this situation both inner peaks vanishes and, in 
addition, the trapping effect disappears. In FIG. [T^ we 
have chosen the initial state as |xj)- For this initial state 
both inner peaks vanishes. However, unlike in FIG. 
the trapping effect is present. Finally, in FIG. we 
display the probability distribution for the initial state 
Ix^)- In such a case the outer peaks vanishes. Moreover, 
the trapping effect is absent, similar to the FIG. [T^ 

Let us now determine the sets of states for which both 
peaks of either the slower walk or the faster walk disap¬ 
pear. To achieve this the weight function has to be of 
the form 



FIG. 11. (Color online) Probability distribution after 100 
steps of the five-state Wigner walk. The initial coin state was 
chosen as IXi")- The coin parameter is p = 0.6 



X 


FIG. 12. (Color online) Probability distribution after 100 
steps of the five-state Wigner walk. The initial coin state was 
chosen as Ixj")- The coin parameter is p = 0.3 



FIG. 13. (Color online) Probability distribution after 100 
steps of the five-state Wigner walk. The initial coin state was 
chosen as \xt)- The coin parameter is p = 0.8 


= (p2 _ ^ ^ ^^2) 
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FIG. 14. (Color online) Probability distribution after 100 
steps of the five-state Wigner walk. The initial coin state was 
chosen as |x^). The coin parameter is p = 0.5 


for some arbitrary a, 5, c. This is satisfied provided that 

= 0 , 


= 0 . 


(34) 


Concerning the slower walk described by the limit density 
from the exp licit form of the weight function 

A4(24)(^;) 

given in ( |C3[ ) we obtain the conditions 
hf +hfho = 0, 

which are satisfied for 


h+ = h^ = 0. 


(35) 


Thus we have a three-dimensional subspace of initial coin 
states 

= holxo) + K\xi) + h^lxt) 

for which both divergencies in disappear. 

For the faster walk described by the limit density 
jy(2,2)(.j^) condition (34) for the weight leads 

us to the relations 

h^{V3ho — h^) + h^{V3ho — h^) =0, 

|y3/io-4|2 + 4|/irP =0. 

This is satisfied provided that 


can eliminate all four divergencies in both densities. The 
state for which this situation occurs is given by 

\^c) = ^\xo) + ^\xt) ■ (37) 

We point out that this is the only coin state for which 
the probability distribution of the five-state Wigner walk 
has only one peak. This peak arises due to the trapping 
effect which we address in the following subsection. We 
illustrate this feature in FIG. fT^ 



X 


FIG. 15. (Color online) Probability distribution after 100 
steps of the five-state Wigner walk. The coin parameter is 
p = 0.5. The initial coin state was chosen according to (371, 
which satisfies both conditions for eliminations of peaks of 
the slower walk (351 and the faster walk (361 simultaneously. 
The resulting probability distribution has only one peak at 
the origin due to the trapping effect. 


Finally, we note that the slower walk described by the 
limit density can vanish completely. Indeed, for 

the state 


\i’c) = l\xo)-^\xt), ( 38 ) 

(2 1 ) 

all terms M] ’ are equal to zero. In such a case the 
spreading of the five-state Wigner walk is described only 
by the limit density We illustrate this effect in 

FiG.na 


B. Trapping probability 


h- =0, h+ = Viho. (36) 

Hence, we find a three-dimensional subspace of initial 
states 

IV’c^) = ^olxo) + K\xt) + V^holx'^) + 


for which the density (v) has no divergencies. 

In contrast to the four-state walk which we have 
treated in the previous Section, it is now possible to sat¬ 
isfy both conditions (35) and (36) simultaneously, i.e. we 


Let us now turn to the trapping probability. We leave 
the details of calculations for the Appendix and fo¬ 
cus on the dependence of the trapping probability on the 
initial coin state. As we discuss in more detail in the Ap- 
pendixj^ the trapping effect does not occur when the ini¬ 
tial state is a linear combination of vectors |x)") and |xl)- 
This feature was illustrated already in FIGs. and [T4| 

It implies that in the suitable basis {|xo)) Ixf)) Ix^} the 

trapping probability depends only on three amplitudes of 
the initial coin state, namely Hq, hi and ■ However, 
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FIG. 16. (Color online) Probability distribution after 100 
steps of the five state Wigner walk. The coin parameter is 
p — 0.7. The initial state was chosen according to (38l. For 
this particular initial coin state the limit density de¬ 

scribing the slower walk vanishes. 


( 2 ) 

the explicit form of (x) is still rather involved. Nev¬ 
ertheless, this can be overcome by an additional change 
of basis in the subspace spanned by {|xo)) IXi")) Ix^} 
which affects the trapping effect. Similarly as for the 
three-state model, which we have treated in Section EYi 
there are coin states for which the trapping effect appears 
only for a; > 0, respectively x < 0. We find that the state 
resulting in trapping of the particle only at non-negative 
positions is given by 


IA+)-y^lxo) + ^lxt) + ^lx})- 

The second state, which traps the particle only at posi¬ 
tions X < 0 is orthogonal to jA’*') and reads 

|A ) = ]/^l^o)-^lxt} + ^IX2)- 

We find that the orthogonal complement of |A^) within 
the subspace spanned by {|xo), Ixf), Ix^} is given by 
the vector (381, which we now denote as |Ao). The triplet 
{|Ao), IA^) {forms an orthonormal basis in the subspace 
affecting the trapping effect, and together with the vec¬ 
tors Ixr) Ix^) if forms an orthonormal basis of the 
whole coin space. When we express the initial state in 
the new basis according to 

ll/'c) = ^o|Ao) + ^^lA"*")-b ^ |A ) +Ixi) + ^2 1X2 )> 
we find that the trapping probability for positive x reads 


p^i(2x) = g2 


and similarly for negative x we find 

„2^ 


P. 


( 2 ) 


(2x) = ^ ( 1^0 + f(x)ri^ + . 


(40) 


Here we have used the notation 


^ ^ + 2|a;|\/l-p2) . 


Directly at the origin the trapping probability has a more 
complicated form 

pS>(0) = ?T(;b!lQ2(i4P + ii_H + 


+ gQ^I^+ + Lp -b 


2 - 




V6(2-p2+ 


8p2 


Q" X 


^ {i^+ + ^-)^0 + (^-1- + L)/o) . 


-Q\lo\^- 


(41) 


The results (n59|) and (40) indicate that the decay of the 


trapping probability with distance from the origin is not 
purely exponential like for the three-state Wigner walk 
(23). Nevertheless, the correction to the exponential de¬ 


cay becomes negligible for large x. The state |Ao) is an ex¬ 
ception, since for 1^ = 0 the terms involving the position- 
dependent function /(x) vanish and the behaviour of the 
trapping probability is exactly exponential. Moreover, 
the decay rate determined by Q is the same as for the 
three-state Wigner walk. 

For illustration, we display in FIG. [m the probability 
distribution of the five-state Wigner walk with the ini¬ 
tial coin state |A+). Clearly, the trapping effect appears 
only for x > 0. Moreover, the plateau formed by the first 
three point indicates that the decay of the trapping prob¬ 
ability deviates from a pure exponential, in accordance 
with (3^. 


VII. CONCLUSIONS 

We have investigated in detail quantum walk models 
on a line with coin operators given by Wigner rotation 
matrices We have shown that this three- 

parameter set of walks is equivalent to a single parameter 
family as far as the shape of the probability dis¬ 

tribution and its moments are concerned. The parameter 
p has a simple physical interpretation since it determines 
the position of the peaks in the probability distribution 
and its width. 

Next, we have simplified the results of [22] by turn¬ 
ing to a more suitable basis of the coin space. Unlike 
for the three-state Grover walk |3^, the suitable basis is 
not directly given by the eigenvectors of the coin oper¬ 
ator, however, we have found a recipe which allowed us 
to construct the suitable basis from them. We presented 
the explicit form of the limit densities for Wigner walks 
up to j = 2. Expressing them in the suitable basis al¬ 
lowed us to identify various interesting regimes which are 
otherwise hidden in the standard basis description. As 
an example, we have shown that the number of peaks 
in the probability distribution can be reduced to one. 


2p4 
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FIG. 17. (Color online) Probability distribution after 10000 
steps of the five-state Wigner walk. The initial state was 
chosen as lA"*"). The coin parameter is p = 0.6. Only a small 
vicinity of the origin is displayed to emphasize that for |A^) 
the trapping effect appears only for x > 0. Higher number 
of steps in comparison with other figures was chosen so that 
the trapping effect is sufficiently pronounced further from the 
origin. Notice the plateau formed by the first three points 
exemplifyin g th e fact that for |A'^) the decay of the trapping 
probability (391 is not purely exponential. 


Initial state 

Properties 

Figure 

Ixo) 

density tends to zero at the origin 

FIG.|2] 

lx+) 

only one peak due to the trapping 

FIG.|^ 

\x~) 

no trapping effect 

FIG.|^ 

|A-) 

no trapping on the positive half-line 

FIG.|| 


TABLE I. Summary of figures for the three-state model. 


Moreover, the suitable basis reveals that p is a scaling 
parameter, since varying the value of p while keeping the 
same amplitudes of the initial state in the suitable basis 
simply changes the width of the probability distribution 
without distorting its shape. 

The Wigner walks with integer j show the trapping 
effect, which was not explicitly evaluated in [32] ■ We 
presented the explicit form of the trapping probability 
for j = 1 and j = 2. Our results show that the trapping 
probability can be highly asymmetric and, moreover, is 
not purely exponential for j = 2. 

For clarity the Tables jT] |ll| and |lTT| summarize the main 
effects and its graphical representation in the text for 
three-, four- and hve-state Wigner walks. 

Based on our explicit results for Wigner walks up to 
j = 2 we make the following conjecture on the suitable 
basis for models with greater j. For half-integer j the 


Initial state 

Properties 

Figure 

Ix^) 

only one peak in the slower walk 
only one peak in the faster walk 

FIG.[^ 

Ix^) 

only one peak in the slower walk 
only one peak in the faster walk 

FIG.[^ 

V3ixr) - ix^) 

no peaks in the slower walk 
only one peak in the faster walk 

FIG.[^ 

Ix^) + ^1X2+) 

no peaks in the faster walk 
only one peak in the slower walk 

FIG.[^ 


TABLE 11. Summary of figures for the four-state model. 


Initial state 

Properties 

Figure 

Ixo) 

no peaks in the slower walk 
density tends to zero at the origin 

FIG. 10 

\xt) 

no peaks in the faster walk 

FIG. 11 

\xi) 

no peaks in the slower walk 
no trapping effect 

FIG. 12 

\xt) 

no peaks in the slower walk 

FIG. 13 

1X2 ) 

no peaks in the faster walk 
no trapping effect 

FIG. 14 

Ixo) + ^1X2+) 

only one peak due to the trapping 

FIG. 15 

Ixo) - ^ 1 X 2 +) 

slower walk vanishes 

FIG. 16 

|A+) 

no trapping on the negative half-line 

FIG. 17 


TABLE III. Summary of figures for the five-state model. 


coin operator has pairs of eigenvectors 

lipn) corresponding to eigenvalues of the form . The 
suitable basis {|Xn)} can be constructed according to the 
recipe 

IXn ) = (^ (e”* " l^/’n) + , 

\Xn) = (e”* " I'^n) - e* = IV'n )) ■ 

For integer j there is an additional eigenvector I'i/^o) = 
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Ixo) with eigenvalue 1, which we include as the last vector 
of the suitable basis. 

The vectors of the suitable basis were selected among 
those which result in non-generic probability distribu¬ 
tion of the Wigner walk. They point out the extremal 
regimes in which these quantum walks can be operated. 
Moreover, they indicate ways in which wave-packet could 
be shaped using quantum walks with higher-dimensional 
coins. We assume that this approach will be useful in 
analyzing the properties of other quantum walk models, 
especially on higher dimensional lattices. 
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Appendix A: Three-state model 
1. Suitable basis 

For the three-state model the explicit form of the eigen¬ 
vectors of the coin (p) reads 


|V^±) = -(±*|l) + v^|0)T*|-l)), 


basis (qi) and the suitable basis (hi) 


|V^o) = ^(|l) + |-l)). 

The eigenvectors satisfy the equations 


(Al) 


qi-^ho + 

qo = + V 

_ 


J 

'l-p2 

V 

2 ‘ 


— p'^h~ 


'l-p^. 

V 

2 ‘ 


h+ - -^h-, 
72 


h'^ H— . 

72 


(A2) 


2. Weight function 

The individual terms of the weight (v) are in the 

standard basis of the coin space given by [42] 

^ (kiP + 2|goP + Ig-ip , 

M^l)=-|gi|2 + |5_i|2 + 


71 - P 

72p 


-{qiqo + qiqo + qoq-i + qoq-i), 


M 


“ 2|(7oP + k-ip) - 


71 — , — _ — ^ , 

——{qiQo + QiQo — QoQ-i — qoQ-i) + 


2-p7 — , - ^ 

- 2p2 (gig-i +9i<?-i)- 


Using the relations (A2) we find that in the suitable basis 
{Ixo)) Ix^)} fhs weight function simplifies into 

= 7+P + 7-7 

P 


= Fo), 


3. Trapping probability 


where the phase tp is given by 

ip = arccos (2p^ — 1). 

The suitable basis is then constructed according to the 
formula ( |2l] ). We find that the vectors of the suitable 
basis {IxoTrfx^)} are in terms of the standard basis given 
by 


|Xo) = ^(|l) + |-l)), 


lx+) = \/^|i) + p|o)- 


1 - p2 


- 1 ), 


2 

lx") = -^|i) + Ti^io) + ^1 -1). 

This leads us to the following relations between the coef¬ 
ficients of the initial state when expressed in the standard 


The trapping effect arises from the fact that the evo¬ 
lution operator of the three-state Wigner walk has an 
infinitely degenerate eigenvalue 1. The overlap of the 
corresponding localized eigenvectors with the initial state 
leads to the non-vanishing trapping probability (10). The 


overlap is simple to evaluate in the Fourier picture m- 
Indeed, the Fourier transformation diagonalizes the step 
operator ([^ and the evolution operator Q reduces into 
a 3 X 3 matrix 


Uik) = Diag(e 


2^k •F^)(p), 


where k is & continuous quasi-momentum ranging from 0 
to 2tt. The matrix U{k) has an fc-independent eigenvalue 
1 with the eigenvector 


v{k) = 


1 


7271 - p 2 

(1 - 6^*'=) p 


74 - 2 p 2 (l + cos2fc) 1 71 


(A3) 
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The overlap of the eigenvector with the Fourier trans¬ 
formed initial state 


tpc = {qi,qo,q-i) 


T 


(A4) 


where the factor Q is given in (24). Finally, the trapping 
probability Px\x) is given by the square norm of the 
amplitude (A5). The final result is presented in (23) 


Appendix B: Four-state model 


where qi are expressed in (A2), yields the limiting prob¬ 
ability amplitude at position x 


2it 

= J (A5) 

0 


Direct evaluation of the scalar product {v{k),ipc) reveals 
that it vanishes when the initial coin state is |x“). Hence, 
for this particular state the trapping effect does not oc¬ 
cur. For other initial states the limiting amplitude can 
be decomposed into integrals of the form 


1. Suitable basis 


For the four-state model the coin operator 
has two pairs of eigenvectors 

\^t) = ;^(^|3/2) T *|l/2) + I - 1/2) T *V3| - 3/2)), 

\^2) = /^(|3/2) ± *V^|l/2) - 731 - 1/2) T *1 - 3/2)), 

(Bl) 


with conjugated pairs of eigenvalues 


ZTT 

X(i)(x) = j 


^—ixk 


47r(2 — _|_ cosfc)) 


dk. 


Using the substitution z = one can turn 
to a contour integral over a unit circle in the complex 
plane, which can be evaluated explicitly by the method 
of residues. We find that the result reads 




Ql"l 

4v/r^’ 


The phases of the eigenvalues are given by 


Lpi = arccos p. 


arccos (p(4p^ — 3)) 


0<p< 2 


P2 = 


27r — arccos (p(4p^ ~ 3)) , | < p < 1 


The suitable basis is then constructed according to (251. 


Using the explicit form of the eigenvectors (Bl) we 


can obtain direct relation between the standard basis 
{|3/2), |l/2), I — 1/2), I — 3/2)} and the suitable basis 
{|x?)) IxJ)}- From this we find the following correspon¬ 
dence between the coefficients of the initial state in the 
standard basis {qi) and in the suitable basis {hi) 


_ 73(1 + P) , 73(1 ~ P) f,- , t/l + P /1 r, \l,+ , 71 ~ P /1 , r, \u- 

^^ 3/2 “ -TTS —^1 -UTS —^1 -„ /t; (1 ~ 2 p)/i 2 H-^ /p; (1 + 2 p)h 2 , 


qi/2 


272 
yl - p 
272 


272 


272 


272 


ht + + 2p)ht - ^^^^^(1 - 2p)h-,, 


272 


272 


272 


+ 2p)/iJ -I- ^ /n^(l ~ 2p)/i2 . 


9-3/2 — 




2^2 


2V2 


2^2 


2V2 


(B2) 


2. Weight function 


The weight functions ^( 3 / 2 , 372 )four-state model in terms of the standard basis 

coefficients qi were given explicitly in |42j , we do not reproduce it since it is rather long. With the help of the relations 
(B2) we can express them in terms of the suitable basis amplitudes hi. We find that the individual terms of the weight 
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function obtain the form 

M(3/2.1/2)^|/^+|2^|/^-|2^ 


Mi' 


(3/24/2) 


M 


M 


(3/24/2) 


(3/24/2) 


1 
P 

_^r 

4p2 




- \Kf) + ^ihth+ + h+ht - - h^h^) 


V3(|/i+p + \hi P - - 1^2 n - + hth+ + /12 ) 


(3|/ii^P - 3\h^ P + P - |/i2 n + Vsihf h+ + / 1 +/ 1 + - hi -h^h^) 


^ 4p3 

Similarly, the coefficients of the weight function for (v) read 


MX 


pi/2,^/2) a/s 


— 2p i^i^2 “b ^t^2 ^1 ^2 ^1 ^2 )i 


M 


M 


(3/2.3/2) 


(3/2.3/2) 


^ r 

4p2 


4p3 


y3(ih+|2 + |hrp - |h+|2 - |h-n - (h+h+ + h+h+ + h-h- + h^h-) 


(3|h+|2 - 3|h-|2 + |h+p - |h-n + V3(h+h+ + h+h+ - h-h- - h-h^) 


(B3) 


(B4) 


Appendix C: Five-state model 


1. Suitable basis 


The eigenvectors of the coin operator M^\p) for the 
five-state model are given by 


IV'2^) = \{\2) ± 2^\l) - 7610) T 2z| - 1) + I - 2)) 

7o) = y||2) + i|0) + y||-2). 


They satisfy the eigenvalue equations 

R^‘^\pMo) = l/’o), 
with phases determined by 
tpi = arccos {2p^ — l) , 

arccos (Sp'^ — 8 p^ + l) ) 0 < p < ;^ 


<P2 = 


27r — arccos ( 8 p^ — 8 p^ + l) 1 ^ ^ ^ ^ 


The suitable basis {|xo), Ixf), lx^)} for the five-state 
model is then constructed according to the formula ([33 1. 


Using the explicit form of the eigenvectors ( jClj ) we can 
obtain the relation between the vectors of the standard 
’ basis and the suitable basis, which leads us to the corre¬ 

spondence between the coefficients of the initial state in 
(Cl) the standard basis and the suitable basis 


® “ 2 ^'*" 
® " 


1 — p2 


h'^ - j=h 


p , l- 2 p\+ , „,/l-p\_ 


-K + p 


2 72 '^' 272 

1 _ 1)2 ,_ 1 _ Pd2 

^ -^r + p72(i - p'^)ht —7^^/’ 


2 ^2 > 


10 


= 4^0 - ^(1 - 2p^)h+ - p73(l-p2)/i2 , 




72 


-hi - Pv/2(1 - p2)hj -b 


9-2 = 


I^-P\+ , P u- , 1-2P%+ , 


272 


hr, — 


-h^ -\ --b 

2 ^ 72 ^ 272 


TiJ -bp 


hrf . 

2 2 


(C2) 
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2. Weight function 


The weight functions and for the five-state model in the standard basis of the coin space can 

be evaluated using the procedure of [42]. Using the relations (C2| we can transform them into the suitable basis. We 
find that the coefficients of the weight function obtain the form 


M 


0^’^^ — + 1^1 


= 


+ |2 , „-,2 

1 r.^— TTT 


P L 
1 r 




= - 3 

p-^ 


3|hoP — — \hi p -|- -|- |/i2 P -I- ■\/3(hoh^ -|- h^h^) 


2(h^^2 /i2 ) “t“ h-^ -t- '\/3(ho/i]^ -t- ) , 

-4\ht\^+ \h+\^+ V3{hoJ4+7^h+) . 

Similarly, the coefficients of the weight function A4^^’^'^{v) in the suitable basis are easily found to be 




,(2.2) 


= \ht\^ + \h 


2 I I 


Mr^ = -- 


=> = 4 






|/ l +|2 + |/ l -|2 _ |/ i +|2 _ |/ i -|2 _ + hoh +) 


2(/i^ /12 “f /z2 ) 4“ ) -f '\/3(ho/i3 -1- hoh-^ ) 

(3|ho|^ - 4|h+|2 + |h+|2) + y3(/io4 + ^/i+) 


(C3) 


(C4) 


3. Trapping probability 


The trapping probability for the five-state Wigner walk 
can be evaluated using the same method described in 
Appendix [A] for the three-state model. The asymptotic 
value of the amplitude at position x is given by 


27r 




dk 


— ixk 


{v{k),'ipc)vik), 


(C5) 


I 

corresponding to the eigenvalue 1. Its explicit form is 
given by 

/ (1 - p2) \ 


;(fc) = 


n{k) 


V 


- (1 - pVi - 

I (p^ cos(2fc) — 2p^ -I- 1 ) 
1 — e’ 
e' 


” 2 *'=) p^l - p2 

.-2.. (1 _ ^2) 


where we have denoted by n{k) the normalization factor 
which reads 


where ipc is the Fourier transformation of the initial state 


i>C = (92, <?i, <?o, 9 - 1 , 9-2) , 


where the coefficients qi are given by (C21. The vector 
v{k) is the eigenstate of the evolution operator in the 
Fourier picture 


C/(fc) = Diag(e'‘*^ 1, ■ R‘'‘^\p), 


n{k) = y -(2 - p {1 + cos(2fc))). 

One can show by direct calculation that if the initial state 
is a linear combination only of vectors |x]”) and |x^) then 
Ipc is orthogonal to v{k). In such a case, the limiting am¬ 
plitude vanishes and the trapping e ffect does not occur. 
Otherwise, the limiting amplitude (C5) can be decom¬ 
posed into integrals of the form 


ZTT 

I^^\x) = J 


^—ixk 


47r(2 — p2(l - 1 - cos k))' 


-dk. 
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which is turned into the contour integral over a unit circle 
in the complex plane and evaluated by the method of 


residues. We find that the result reads 


16(1-p2)i 


where Q is given by (^24|. The trapping probability 

( 2 ) ' -' 

Px\x) i s th en obtained by the square norm of the am¬ 
plitude (C5|. The final result for positive x, negative x 
and the origin is presented in (39), (40) and (41). 
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